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^ ■ Abstract 
(N 

\Q ' In principle, leading neutrons produced in photoproduction and deep-inelastic scat- 

tering at HERA have the potential to determine the pion structure function, the neutron 
■ absorptive cross section and the form of the pion flux. To explore this potential we com- 

pare theoretical predictions for the and pt spectra of leading neutrons, and the Q 2 
dependence of the cross section, with the existing ZEUS data. 
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1 Introduction 

The original motivation for measuring leading neutrons at HERA was to determine the pion 
structure function F£(x, Q 2 ), assuming that the process jp — > Xn is dominated by 7r-exchange 
and that the 7r-flux is known from the analysis of hadronic data in the triple-Regge region. 
The problem, however, is that the original 7r-flux, given by the simple 7r-exchange diagram, 
is modified by "soft" rescattering effects which are different in deep inelastic scattering and 
hadron-hadron collisions. In hadron-hadron collisions the leading neutron has a much larger 
probability of a secondary interaction, which changes the longitudinal fraction, x^, of the 
proton's momentum carried by the neutron, as well as the neutron's transverse momentum, p t . 
On the contrary, for deep-inelastic scattering the probability of secondary interactions becomes 
negligible; and so knowing the 7r-flux, we can directly extract F£{x, Q 2 ) at large Q 2 . 

At HERA we have the possibility of varying Q 2 . So the increasing role of rescattering can 
be studied as we go to lower Q 2 , and finally to photoproduction. It was shown in Ref. [T], 
that for xl > 0.6, rescattering may just be considered as absorption of the fast neutron. After 
rescattering the neutron migrates to lower xl, but, at large xl, the population due to migration 
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is small and here rescattering acts as absorption. Thus, for large xl, absorptive effects may be 
described in terms of the survival factor, S 2 , of the rapidity gap associated with 7r-exchange. 
S 2 may be calculated in the usual way (see, for example, [2]); its value depends on the inelastic 
neutron cross section. Note that a 7*-initiated process proceeds through a 7* — > qq transition, 
followed by the interaction of the gg-pair with the target. The probability of rescattering 
is therefore determined by the cross section of this qq interaction. The size of the gg-pair 
decreases with increasing Q 2 , and the effect of rescattering becomes negligible at large Q 2 . 
Therefore leading neutron data for a range of Q 2 gives information about the gg-neutron cross 
section, or <r a b s - The preliminary analysis of the data indicates that both F^ and cr a b s are in 
reasonable agreement with theoretical expectations based on the additive quark model. 

The next step is the possibility of studying the precise form of the 7r-flux by measuring the 
p t dependence of the neutron yields. Thus, thanks to the additional parameter Q 2 , it seems 
possible to extract three types of information from one set of leading neutron deep-inelastic 
data: 

(a) the pion structure function F£(x, Q 2 ), 

(b) the (qq) — N absorptive cross section cr a b s , where N = n or p, 

(c) the form of the 7r-flux; the normalization of the 7r-flux is given by the known G-nNN 
coupling constant, when the pion is close to its mass shell. 

Unfortunately, at present, the p t dependence measured at HERA is flatter than that expected 
from theoretical models — that is from any reasonable parametrization of the pion form factor. 

In Section 2 we compare our cross section predictions for leading neutrons, as a function of 
xl, with the HERA data, and in Section 3 we study the p t dependence. In order to account 
for the flatter p t dependence of the data, in Section 4, we study the possible role of p and a 2 
exchange contributions. We consider both photoproduction and deep-inelastic production of 
leading neutrons. 



2 Photo- and DIS-production as a function of xl 

If we assume 7r-exchange dominance, then the inclusive production of leading neutrons, jp — ► 
Xn of Fig. I^a), is given by the triple- Regge diagram shown in Fig. db). We have 

da(^p ->■ Xn) 2 G^+pn (~t) 2 \l-2a*(t) _tot/ n/f 2\ /-,\ 

dx^dt " S Tg^~ (t - ml) 2 F W (1 " XL) ^ (M) ' (1) 

where the coefficient of is called the pion flux. The pion trajectory, a n (t) = a' n (t — m 2 ), is 
taken to have slope a' n ~ 1 GeV -2 , and the 7r + pn coupling constant is G 2 +pn /8n = 13.75 0. 
The invariant mass M of the produced system X is given by M 2 ~ s(l — Xl). F(t) is the form 
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Figure 1: (a) The pion-exchange amplitude and (b) the corresponding dominant triple- Regge 
contribution to the cross section of the inclusive production of leading neutrons at HERA, 
7P — > Xn. 



factor resulting from the pion-nucleon and mrP vertices with off-mass-shell pions, see Fig.^b). 
The survival factor S 2 , which takes into account absorptive corrections, depends on x^ and p t 
of the leading neutron. The calculation of S 2 is outlined in the Appendix. 

The cross section of the 77r- interaction, 0"*°*, and the pion structure function, F£, are the 
quantities measured in photoproduction and deep inelastic scattering respectively, where 

47T 2 « 

~Q 



<t = -TyT F 2- (2) 



We use the additive quark model to obtain theoretical estimates, assuming for photoproduction 
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and for deep-inelastic scattering 1 



f^x,q 2 ) = \f p 2 {\x,q 2 ). (4) 

We rescale the Bjorken x in order to have the same energy for the 7*-valence q interaction. 
Another possibility, which we will discuss, is to simultaneously rescale Q by the ratio of the 
pion and proton radii. It was shown in Ref. pQ, that if we take a reasonable value of the 
neutron absorption cross section 2 then this approach satisfactorily describes the ZEUS data 
for the photoproduction of leading neutrons at large x^. The description, updated for the new 
experimental cuts used in 4\, is shown in Fig. |21 From the figure we see that the absorptive 
corrections reduce the cross section, given simply by Reggeised pion exchange, by a factor S 2 , 
averaged over p 2 , of about 0.5 independent of xl ■ 

1 Unfortunately, the present parametrizations of the parton distributions of the pion are unreliable in the low 
x region of interest. Therefore we take (@J. 

2 The value taken was motivated by the p-dominance model of the photon. 
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Figure 2: The predictions for the Xl spectra of photoproduced leading neutrons compared 
with preliminary ZEUS data [4J; only the systematic errors on the data points are indicated, 
as these dominate the statistical errors. The dotted, dashed and lower continuous curves are 
respectively the results assuming first only Reggeised n exchange, then including absorptive 
effects, and finally allowing for migration; the calculation is described in [T], updated here to 
allow for the different experimental cuts. The upper continuous curve corresponds to including 
p and a,2 exchange contributions, as well as 7r exchange, as described in Section 4. 
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From the theoretical point of view, it would be best to observe leading neutrons produced 
in DIS at very large Q 2 where the rescattering absorptive corrections are negligible; and to 
measure F^ in a most direct and clear way. Unfortunately, the event rate at large Q 2 is limited. 
The ZEUS preliminary data jl] correspond to Q 2 > 2 GeV 2 , with an average, (Q 2 ), of 16 GeV 2 , 
so we cannot neglect absorption even in the DIS data sample. To be precise we have to integrate 
over the size of the qq pair produced by the photon, starting from the hadronic/confmement 
scale ~ Aq CD ~ 0.1 GeV 2 up to (Q 2 ). The part of the cross section originating from a small 
size qq pair will have negligible absorption, while a large size pair will be strongly affected by 
rescattering and will suffer an S 2 suppression. This prescription was implemented explicitly in 
Refs. (HI U\ ■ Here we adopt a simplified approach assuming that the part of F% measured at 
the initial hadron scale 3 enters with the same absorptive factor as in photoproduction, while 
the remainder of F£, which is generated by DGLAP evolution with strong fc t -ordering, that is 
by a small size qq pair, is taken to have S 2 = 1. 

Since the absorptive effects in the present DIS data are not negligible, the ratio of the 
probabilities of observing leading neutrons in photoproduction compared to DIS production 
will be closer to 1 than the ratio of the curves shown with and without absorption in Fig. El 
The ratio is plotted in Fig. El for two different values of Q 2 . The predictions are calculated with 
the same cuts as the ZEUS leading neutron data jS] . The growth of the ratio as xl — > 1 reflects 
the fact that the energy (that is, the x) dependence of DIS production, a ~ x~ x , is much 
steeper than for photoproduction. Since x in the interaction is smaller by a factor (1 — xl) 
than in the proton structure function used to normalise the neutron yields, (da / dx l) / cr , this 
gives an extra factor 

(1_ Xl )A(o)-a(qVo) (5) 

in the ratio plotted 4 in Fig. El For xl ~ 0.7, where this factor is not so important, the ratio 
R decreases with increasing Q 2 since at larger Q 2 we have smaller absorption. Of course the 
prediction of the ratio depends on the value chosen for cr a b s and the ansatz used for F£. For 
example, if we used the rescaled value of Q 2 in eq. (|IJ) then we would increase the predicted 
ratio by 10 — 20%. Turning the argument around, we see that if the ratio is precisely measured 
for different Q 2 values, then it would be possible to determine both F% and cr a b s - 

3 The p t dependence of leading neutrons 

The predicted p t dependence of leading neutrons is shown in Fig. for different values of x^. 
We see that the distributions do not have exactly exponential form, as may be anticipated from 
(JTj). However the departure from this form is not large, and the experimental data are usually 
discussed in terms of an average slope b, where da/dp 2 ~ exp(— bp 2 ). Note that, in spite of the 
factor of t in the numerator of (QJ), the distributions do not vanish at p t = 0, due to 

Uin - - (1 - x L ) 2 m 2 N /x L ; (6) 

3 We take the scale to be m 2 = 0.6 GeV 2 . 
4 We use MRST2001 LO partons 0. 
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Figure 3: The ratio, R = photoprod./DIS, of leading neutrons from photoproduction, (da(Q 2 ~ 
0)/cIxl)/o'(0), and DIS production, (da(Q 2 ) / dxi) / o{Q 2 ) as a function of xl for two different 
intervals of Q 2 , namely 4 < Q 2 < 7 GeV 2 and Q 2 > 13 GeV 2 , compared with predictions for 
Q 2 = 5.3 GeV 2 (upper curve for xl < 0.9) and Q 2 = 51 A GeV 2 (lower curve). The predictions 
include the p and 02 exchange contributions, as well as 7r-exchange, and the rescaled value of 
Q 2 (eqs. (JBJ), Q), as described in Section 4. The points are preliminary ZEUS data [3]. 
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where t = — (p 2 /^) +£ m in- The calculated distributions qualitatively agree with the data. The 
slope, b, grows with xl, as can been seen in Fig. |5J This arises from the factor (1 — x^) -20 ^ 
in ((H), and reflects the Reggeization of the pion. The predicted slope for photoproduction of 
leading neutrons is a bit larger than for DIS production, see Fig. |3 since strong absorption at 
small impact parameters, pt, effectively pushes the distribution into the larger pj- region. 

Unfortunately, Fig. shows that the predicted p t distribution is too steep in comparison 
with the data. This is a general property of any model, see Of course, there is a possibility 
to modify eq. (JTJ), using a more complicated signature factor. If we replace the pion propagator, 
l/(m 2 r — t), by the usual Reggeon signature factor, 7ra^/(2 sin(7ro; 7r (t)/2)), then the slope is 
a bit smaller, but only by 0.2 — 0.4 GeV -2 . In some studies, motivated by Veneziano-type 
models, the signature factor is written in terms of the T function. Such a function, r(— a^it)), 
grows at large negative t, but, first, it does not help in the region of interest, —t < 0.5 GeV 2 , 
and, second, this growth reflects the contribution of heavier mesons which lie on daughter 
trajectories. The simplest possibility, to reconcile the calculated slopes with the measure values, 
is to take a negative slope of about 1 GeV -2 in the form factor F(t) in (fT|). but this is physically 
unacceptable. 



4 Including p and a<i exchange contributions 

This persistent discrepancy between the measured and predicted p r slopes of leading neutrons 
from 7r-exchange leads us to investigate the contributions of secondary (p, a 2 ) Reggeon ex- 
changes. Note that the spin-flip a 2 -exchange amplitude interfers with the 7r-exchange ampli- 
tude. These contributions were discussed in Ref. JS], where it was stated that they never exceed 
10% of the 7r-exchange induced cross section. But, there, the authors concentrated on very low 
p t , to be close to the 7r pole, and did not allow for the spin-flip p and a 2 nucleon- Reggeon 
vertex. The exchange of the p-trajectory (but not of the a 2 -trajectory) was considered also in 
|10j . There the p contribution was estimated to be about 20% of the whole cross section for 
events with a leading neutron with xl = 0.7 — 0.8. As a result the authors pointed out that 
including p exchange will lead to a smaller slope in the p t distribution of leading neutrons. 

For p-exchange, spin flip is very large [TT] , 



V p &ip g Pt 

^non-flip 2m N 



(7) 



To evaluate the effect of p, a 2 exchange, we use the parametrization of Ref. ^1] for the p meson, 
assume p — a 2 exchange degeneracy, and assume the same additive quark model relations for 
the 7p and 7a 2 cross sections, which we used for 771- , that is (jHJ) and (jlj). Clearly these extra 
exchange contributions will enlarge the cross section for leading neutron production, mainly at 
the larger p t values, and less as p t — > 0, which will lead to a smaller average slope. Indeed the 
calculated slope is in more agreement with the data, but now the probability to observe leading 
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Figure 4: The p 2 distributions of leading neutrons for four different values of x^. The lower 
two curves correspond to photoproduction (PHP) and production by deep inelastic scattering 
(DIS) with Q 2 = 16 GeV 2 . The upper curve is obtained assuming that there is no absorption. 
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Figure 5: The slope, b, of the pi distribution of leading neutrons produced in DIS as a function of 
xl- The data are taken from [I]. The upper curve corresponds to ir exchange, with absorption 
and migration effects included as in Ref. pQ, but with the same kinematic cuts as used to 
obtain the data jl]. The lower curve is calculated including p and a 2 exchanges, in addition 
to 7r exchange, as described in Section 4. The decrease of the slope for xl ^ 0.85 reflects the 
vanishing of t m i n as xl — > 1. 
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Figure 6: The difference in the slopes of the pt-distributions of leading neutrons coming from 
photoproduction and from deep inelastic scattering. We compare the preliminary ZEUS data 
jl] with the calculations based on n exchange (upper curve) and then including also p and a 2 
exchange (lower curve). 
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neutrons becomes too large. That is by using the same values of the parameters that we took 
in [T], and the additive quark model relations of (jHJ) and (J1J), the calculated cross section is 
nearly 40% above the data. 

So we have better agreement with the data for the slope b, at the expense of now failing to 
describe the size of the cross section. Can we do better? Of course, there are many parameters 
which are not precisely known. In an attempt to achieve a simultaneous description of the 
cross section and the slope of the p\ distribution of leading neutron production, we change 
the parameters within the limits of acceptability. We take the slope of the pion trajectory to 
be a' w = 0.9 rather than 1 GeV -2 . We use the normal Reggeon signature factor of the pion 
7ra4/(2 sin(7ra 7r (t)/2)). We increase the absorption by setting C = 1.6 and <7 tot (Trp) = 34 mb, 
rather than 5 C = 1.3 and cr tot (7rp) = 31 mb as in pQ. This has the effect of reducing 6 S 2 from 
0.5 to 0.4. In addition we should be more careful with the form of the additive quark model 
relations, (jHJ) and (@J), for the 7-meson cross sections. Let us discuss the perturbative QCD 
expectations. At lowest order in as, the cross sections are proportional to the radius squared 
of the hadron, o ~ a 2 s r 2 [HJ EH E] . If the two colliding particles have quite different radii, 
r a "C r&, the cross section is controlled by the smallest radius. Thus in DIS, the 7-meson cross 
section a ~ 1/Q 2 , as in (j2J). However the interval of evolution of the structure function F 2 
starts from the largest radius r b . Thus we should rescale the value of Q 2 in (JU), and take 

Fr s ° n (x, Q 2 ) = ? Ff Qx, ^Q 2 ) . (8) 

Moreover, we assume that the quark wave functions in the p and 7r mesons are the same (as 
would follow from SU(6) symmetry). We take r 2 ^ = 0.44 fm 2 ^H] and r 2 = 0.76 fm 2 ^H]- On 
the other hand, for photoproduction, where the two radii are comparable, it is reasonable to 
replace r 2 by r m r p . Then we have 

rr tot - — rr tot (Q) 

7— meson 9 TP ' v / 

in the place of ©• 

The predictions obtained with these modifications are shown in Fig.|21and Fig.EJ Indeed we 
do achieve a more satisfactory simultaneous description of the cross section and the p 2 -slopes 
as functions of xl. The results of Fig. H]now become those shown in Fig. The lower pair of 
curves correspond to just the 7r-exchange contribution. By comparing the 7r component with 
the full contribution, we see that, for very low p t , 7r-exchange provides more than 70% of the 
cross section for xl > 0.7, whereas at p 2 ~ 0.3 GeV 2 the (p, (^-exchange contribution starts 
to dominate. Finally, the curves plotted in Fig. El correspond to the modified description. 

5 There is evidence that the factor C, which accounts for the diffractive excitation, should be larger • At 
first sight we might expect that stronger absorption would be generated by including the contribution from the 
"enhanced" diagrams, which arise from rescattering from intermediate partons. However then there would be 
a strong energy dependence for the probability to produce leading neutrons. This is not observed in the data. 
Further discussion is given in Ref. pp. 

6 This is more in line with an earlier calculation which gave a rapidity gap survival factor S* 2 of 0.34 for 
the resolved part of the photon wave function. 
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Figure 7: The p\ distributions of leading neutrons for four different values of xl as in Fig. 
but now including the p and ci2 exchange contributions, with increased absorption, as described 
in the text. For completeness, we also show, by the two lower curves, the individual tt ex- 
change contribution for DIS (dotted) and photoproduction (continuous) using the modified 
prescription. 
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5 Conclusions 



Leading neutron spectra at HERA have a rich structure. The data offer a serious challenge; 
it is difficult to describe the dependence of the cross section on xl and Q 2 , while at the same 
time satisfactorily reproducing the observed p t behaviour. First pure Reggeized 7r-exchange 
fails on all counts. The cross section data indicate the important role of absorptive corrections, 
arising from rescattering. This effect suppresses the yield of leading neutrons by a factor of 
about 0.5 in photoproduction, but less in DIS production. For the same reason the slope of the 
p t distribution in photoproduction is a little larger that in DIS production. The calculation, 
based on Ref. pQ, but updated to account for the new experimental cuts gives a satisfactory 
description of the observed cross section for xl ^ 0.6. 

However, although the difference in slopes for photoproduction and DIS production is rea- 
sonable, the predicted p t distribution for photoproduction falls off more quickly than the data. 
This is a long standing general problem for 7r-exchange models. Indeed, the relatively low 
experimental value of the slope b indicates the presence of secondary (p, a 2 ) Regge spin-flip 
contributions. It is desirable to measure the slope in different intervals of pt since it offers the 
chance to separate 7r from p and exchanges. Recall that for it exchange, the vanishing of the 
amplitude as t — > is almost compensated by the n propagator l/(t — ml). On the other hand 
spin-flip p-exchange is indeed proportional to t. Therefore in the largest xl interval, where t m i n 
becomes small enough the presence of a significant p contribution would produce a "kinematic 
turnover" in the slope. In fact, neglecting t m i n , the spin-flip p contribution leads to a negative 
slope as pt — ► 0. This behaviour explains the different behaviour of the two curves shown in 
Fig. El as xl — > 1. 

Although the presence of the secondary trajectories improve the description of the abso- 
lute slope, they disturb the satisfactory description of the cross section obtained by just n 
exchange. We were required to look in more detail at the calculation of the survival factor. A 
physically-motivated recalculation of absorption yielded S 2 = 0.4. In this way we achieved a 
more satisfactory overall description of the leading neutron spectra. The original and the mod- 
ified pion structure functions, F£(x, Q 2 ), that is (JH) and (JHJ), are plotted in Fig. |H]as functions 
of Q 2 for three different x values, together with predictions obtained from parton distributions 
of the pion. 

Already an initial analysis of the preliminary ZEUS data have yielded interesting results. 
To obtain a reasonable description of the all features of the leading neutron data we are lead 
to select the pion structure function given by the continuous curves in Fig. |H] in preference to 
those given by the dashed curves. Moreover, the pion structure function preferred by the data 
is found to be larger than that obtained using the parton distributions of the pion, see, for 
example, (201 • Here, we should note that the parton distributions of the pion were obtained 
from an analysis of high statistics data, from tt^N experiments for both Drell-Yan and prompt 
photon production, which did not extend below x n ~ 0.2. So extrapolation is necessary to 
obtain the SMRS curves in Fig. |H1 
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In summary, the physically- motivated modification, (|SJ), while not perfectly reproducing 
every feature of the leading neutron data, is clearly a step in the right direction. The analysis 
described here should be regarded as a broad-brush exploratory study. In our additive quark 
model type of approach we have forced the structure function of the pion to mimic that of the 
proton. There is no reason why it should, and so some mismatch with the data is to be expected. 
Indeed we can expect a more detailed analysis of more precise data with a freely parametrised 
pion structure function, F£(x, Q 2 ), will give a quantitative measure of the behaviour of the 
function, as well as the properties of the rapidity gap survival factor, S 2 , that is, of absorptive 
effects. 
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Appendix 

It is simplest to calculate the absorptive corrections in impact parameter, p?, space. For 
convenience, we adopt the notation used in Ref. [I]. We take the Fourier transform of the 
original 'bare' amplitude, A, which corresponds to the cross section given in (0) without the 
S 2 factor. It depends on the impact parameter p nn . The subscript n denotes the 'transverse' 
position of the 7r — 7r— Pomeron vertex in Fig. ^b). That is, p vn is Fourier conjugate of the 
neutron transverse momentum in the bare 7r-exchange amplitude. Then the transverse distance 
between the leading baryon and the incoming photon (or gg-pair) is 

PT = Pyn = |Arn + Pynli (10) 

where p 7?r is the impact parameter for the amplitude describing the interaction of the incoming 
photon with the 'effective' pion (exchanged in the t-channel). For a fixed value of p 77r , the 
suppression of amplitude A is given by the eikonal exp(— fi/2), with 

V(pt) = g| exp(-p 2 T /AB) , (11) 

where the slope B of the p7r amplitude is taken to be 5 GeV -2 . Thus to determine p t we must 
multiply A(p nn ) by the eikonal factor exp(— Q/2), and perform the inverse Fourier transform 

A{x L ,vuM = /^e**** A{x L ,p m )e- a ^f\ (12) 
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Figure 8: The pion structure function, F£{x, Q 2 ), as a function of Q 2 for three different values 
of x. We also show the structure function calculated from the partons of [2U|. On this plot x 
denotes the Bjorken-x of the pion, x n , that is x n — x/(l — Xl). 
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Finally, we integrate over p 77r to obtain the cross section 

da _ f 2 \A(x L ,p t ; p 77r )| 



2 



d 2 p, n ' W ! e" p - /4B , (13) 



dx L dp 2 t J ryn 4nB 



suppressed by rescattering effects. 
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